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Abstract 

We present the complete threshold enhanced predictions in QCD for the total cross section of 
gluino pair production at hadron colliders at next-to-next-to-leading order. Thanks to the compu- 
tation of the required one-loop hard matching coefficients our results are accurate to the next-to- 
next-to-leading logarithm. In a brief phenomenological study we provide predictions for the total 
hadronic cross sections at the LHC and we discuss the uncertainties arising from scale variations 
and the parton distribution functions. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) is an attractive extension |lj[2j of the 
very successful Standard Model (SM) of particle physics. One property of the MSSM is its rich 
spectrum of heavy particles which might be discovered at the LHC if they are lighter than « 2 TeV. 
Hadron colliders are especially appropriate to study color-charged particles. In the MSSM, the 
superpartners of the gluon and the quarks are the gluino, which is a Majorana fermion, and the 
scalar quarks (squarks), respectively. The associated superpartners of the SM particles have the 
same weak isospin, hypercharge and color charge. 

Searches for supersymmetry have been performed at the Tevatron and the LHC with center- 
of-mass (cms) energies of 7 and 8 TeV. Due to its larger energy the best bounds on the masses 
of these hypothetical particles come now from the LHC. The production cross section for gluino 
pairs at the LHC is sizable. It is driven by the large gluon luminosity and it is further enhanced 
due to the large color charge of gluons and gluinos. However, no superpartners of the SM 
particles have been discovered so far. 

A special feature of particle spectra in the constrained MSSM (CMSSM) are the large mass 
differences between squarks, gauginos, and sleptons, allowing cascading decays of the SUSY par- 
ticles. At the Atlas experiment, searches for squarks and gluinos in a CMSSM framework are 
performed by looking for final states with a large number of jets and missing transverse momen- 
tum |3j, additional same sign leptons |4j or b-jets (5). Gluino masses smaller than 840 GeV |3j, 
550 - 700 GeV g), and 600 - 900 GeV (5) are excluded and similar results (6}jl0} are reported by 
the CMS experiment. The LHC bounds discussed above do not apply if the particle spectrum is 
compressed. In such scenarios, if the gluino is mass degenerate with the lightest supersymmetric 
particle and the squarks are decoupled, a lower mass bound of the gluino mass of 500 GeV holds 
(see Ref. pTj for a detailed discussion on how the LHC bounds change). 

Theoretical predictions for gluino pair production up to next-to-leading-order (NLO) in QCD 
have been obtained in Ref. jl2| . The hadronic leading order (LO) and NLO cross sections can be 
evaluated numerically using the program Prospino |13) . As an improvement beyond NLO, the 
threshold enhanced logarithms have been resummed to next-to-leading-logarithmic (NLL) accu- 
racy |!4"fll5| , implying corrections of about 2 — 35% in comparison to the NLO cross section which 
depend on the gluino mass and the chosen parton distributions (PDFs). In Ref. fll6| [, threshold 
effects at NLO in QCD due to remnants of the IS resonance of gluino bound states are discussed 
leading to an enhancement of the complete NLO threshold cross section of 7 — 9% compared to 
the fixed order predictions. Recently, the combined NLL resummation of threshold logarithms 
and the Coulomb corrections for gluino pair production has been studied in Ref. (17) and phe- 
nomenological predictions for cross sections at the LHC have been summarized in Ref. |18|. The 
inclusive cross section for squark-antisquark pair production has been subject to similar improve- 
ments in the past 1 14 19 1 and presently, the corresponding predictions beyond NLO are exact to 



next-to-next-to-leading-logarithmic (NNLL) accuracy [20) . In contrast, the available results for 
gluino pair production are still limited to NLL accuracy, only. 

In this article, we improve the available QCD predictions for gluino pair production to NNLL 
accuracy, putting it on par with the case of squark-antisquark pair production. To that end, we 
compute the missing hard matching coefficients at NLO near threshold. With our new results, 
we are able to provide QCD predictions for the total hadronic cross sections at approximately 
next-to-next-to-leading order (NNLO). These corrections lead to a further increase of the cross 
section of the order of 10% in comparison to the NLO results. As all searches for SUSY particles 
so far have resulted in exclusion limits only, a precise knowledge of the gluino pair production 
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Figure 1: Feynman diagrams for the production of a gluino pair gg via qq annihilation (diagrams (a) - (c)) 
and gg annihilation (diagrams (d) - (f)) at LO. 

cross section in the threshold region is of special interest, because the size of the expected rates 
has a direct impact on the excluded mass range for gluinos. 

The article is organized as follows. In Sec. |2j we recall the basic ingredients of the hadronic 
and partonic production cross sections. In Sec. |3j we review the formalism of threshold resum- 
mation in Mellin space and then proceed to extract the color-decomposed NLO cross section at 
the threshold from known results for gluino-bound state production given in Ref. [16 1. Verify- 
ing a general result of Ref. [21 1, we then calculate the color-decomposed NNLO cross section 
in the threshold limit. Finally, we resum the cross section to NNLL accuracy, matched onto the 
approximated NNLO result. We check our analytic formulas by extracting the color-summed 
one-loop matching constants from Prospino via an appropriate fit in the threshold region. The 
hadronic production cross section is discussed in Sec. [4] The appendices contain useful analyt- 
ical expressions for certain scalar n-point integrals and the expansion coefficients of the general 
resummation formula. 



2 Theoretical Setup 

We study the hadro-production of gluino pairs at the LHC (i.e. the reaction pp — > gg) with its 
partonic sub processes 

gg -> gg, (1) 
W gg, 1 = d,u,s,c,b , (2) 

with the relevant LO Feynman diagrams shown in Fig. [T] 

We focus on the inclusive hadronic cross section of hadro-production of gluino pairs, Cpp^-ggx* 
which is a function of the hadronic cms energy i/s, the gluino mass nig, the squark mass 
(assuming mass degeneracy among the squark flavors), and the renormalization and factorization 
scales, y. r and fit, respectively. In the standard factorization approach of perturbative QCD, it 
reads 
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where the parton luminosities L,y are given as convolutions of the PDFs fi/„ defined through 



Lij(s,s, jij) 



dz 

2 



fi/p (^/'-)///p \ ff> 



(4) 



Here, s denotes the partonic cms energy. As will be discussed below, the treatment of soft-gluon 
exchange in the final-state gluino pair requires the knowledge of the color-decomposed partonic 
cross sections c^-j. Setting ]ij = ji r = \i, the latter is commonly expressed by dimensionless 

scaling functions fjffi in a power series in the strong coupling constant a. s = tx s (}i), 
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r(10) , r(ll)i 



2/r(20) . ,(21), 



f(22) r 2 
'ij 



, (5) 



with Lp t = ln(^ 2 / m|). We use the capital index I to label the admissible SU(3) co i or representations 
of the scattering reactions ^ and Q. The decomposition of the gluon-fusion channel ([TJ into 
irreducible color representations is given by 



8x8 



8 S + 8„ + 10 + 10 + 27 s , 



(6) 



and a suitable basis in terms of the generators of the adjoint representation of the SU(3) co i or can 
be found in [22] (see also |14|[). Likewise, for the quark-antiquark channel we use the color 
basis 



3x3 = l s + 8 s + 8 fl . (7) 
The partonic cross sections in Eq. ^ are recovered after summation over all color structures, 
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(8) 



and, similarly, for the scaling functions in Eq. (pi. The (color-summed) scaling functions at LO 
are given by |12| 
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with the abbreviations 
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(11) 
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The threshold expansion of the LO hard function for gluon-fusion depends only on the dimen- 
sionless variable /3, which is zero at the threshold s = 4m~. The expanded LO function for 
quark-antiquark annihilation depends on /3 and on the ratio r of the squared squark to gluino 
masses. For r = 1, the linear term vanishes. Therefore, small mass differences between the gluino 
and the squark mass lead to a suppression of the qq channel. The gq channel on the other hand 
is absent at tree level. Its NLO contribution at threshold is of the order /3 3 ln(/3) and thus strongly 
suppressed compared to the gg and qq channels. We will therefore ignore its contribution to 
gluino production in the following discussion. However, we include its NLO contribution to the 
total hadronic cross section. 

For the color decomposition of Eq. ^ we find in agreement with Ref. fll4| the following Born 
scaling functions for the reaction dlj, 
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(13) 



(14) 



with Li given in Eq. ([TTJ. Note, that in the threshold limit, only symmetric color representations 
in Eq. p2| contribute, for which we define the normalization factor, 
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for I = {1,8,27}. 



(15) 



For qq annihilation in Eq. ([2]), we obtain the LO scaling functions in the color decomposition 
of Eq. as 
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with L2 given in Eq. pi) . Note that at threshold, only the antisymmetric octet representation of 
the qq channel contributes if the gluino and the squark masses are different. If the gluino and the 
squarks have equal masses, the antisymmetric octet scaling function is vanishing up to C(/3 4 ), 
see Eq. l [T8} , and the symmetric singlet and octet scaling function contribute with the ratio 2 : 5 
at the production threshold. If r ^ 1, the gluino pairs are produced in an S-wave, otherwise in a 
P-wave in that channel. 
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3 Higher order partonic cross sections at the threshold 

At higher orders in QCD, the cross sections develop large threshold logarithms of the type ln(/S) 
in the region s ~ Am 2 , which can be resummed systematically to all orders in perturbation theory. 



Here, we make use of techniques described in 1 23 -291. The resummation is performed in Mellin 



space after introducing moments N with respect to the variable p = Am-/ s of momentum space, 



&(N,ml) = J dpp N - l a{s,m\). 



(19) 



As the threshold limit f> — > corresponds to N — > oo, all terms proportional to powers of 1 / N 
will be discarded. The general resummation formula reads 



% i (N, m 2 ~) = &l , (N, m 2 ~) g\ l (N, m§) exp \c ih i(N + 1)1 + ©(AT 1 In" N) , 



(20) 



where we have suppressed all dependence on the renormalization and factorization scale, }i r 
and fijr. The subscripts ij denote the production channel, where we consider ij = gg,qq. The 

exponent Gy; i contains all large threshold logarithms ln^ N in Mellin-N space, and the resummed 
cross section, as indicated in Eq. ( |20| , is accurate up to terms which vanish as a power for large 
Mellin-N. To NNLL accuracy, G !; - ; i is commonly expanded as 



G !7/I (N) = In(N) -gUX) +gf irl W+a sg U\) 



(21) 



where A = a s fio InN and we abbreviate a s = oc s / (4zr). The functions g 1 -- 1 are derived from the 
double integral over a set of anomalous dimensions (see e.g., 1 29 -31 [), 



G,;,l(N) 



7 N-l 



dz 



4m2(l-z) 2 



dq 2 



Aj(a s (q 2 )) +Aj(a s (q 2 



(22) 



+D i; - I ( as (4 m |(l-z) 2 ))|. 

Here, the cusp anomalous-dimension A\ refers to initial-state collinear gluon radiation, while any 
large-angle soft gluon radiation is contained in the function D,^ \, which splits into the functions 



Di(a s )+ D ; (« s )j + D m {oc s ), 



(23) 



for initial- and final-state radiation, where D, can be taken from threshold resummation for the 
Drell-Yan process or for Higgs production in gluon-fusion. The perturbative expansions for the 
anomalous dimensions reads 



Ai(cis) 



my* 







(24) 



(same for D, (a s ) etc.) and the expansion coefficients A*p and are both known to third order in 



a s from Refs. |32 



the final state depends on 



33|and|34 



35 



, respectively. The function D.) — j due to soft gluon emission in 



the SU(3) co i or representation of the final-state gluino pair and results up 
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to second order in a s are given in Ref. |3l) for heavy final states in arbitrary color representations^ 
(see also 1361). This suffices to compute the functions g\j\ in Eq. (Ell to NNLL accuracy even 



with the dependence on the }i r and ji^ separated (for the computation see, e.g., Refs. |29 37|). The 
explicit expression for g^ 1 can be read off from Eq. (A.5) of Ref. |30|, for g^ 1 from Eq. (A. 7) and 
for g|j from Eq. (A.9) of that reference with the replacements A^ 1 '-^ A^/ff, dJ /} n(0 ' a{l) 
D 



d1^j//3q' (note the sign convention for Dqq in |]30|), and f>\ 



Pi/p>l +1 \ where 



QQ ~^ "gg.i'ro 

denote the well-known QCD beta-function coefficients in the normalization ( |24) |. 

As a last remaining step in achieving resummed predictions to NNLL accuracy in QCD, one 



has to extract the process-dependent matching constants g° j in Eq. (20 1. These consist of the 



hard coefficients gj?. i(« s ) multiplied by Coulomb coefficients g®j I (« S ,N), which also account for 
the interference of Coulomb exchange with hard contributions and soft radiation. A perturbative 
expansion in analogy to Eq. d24) yields 



,(0),C, 
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(25) 
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(N)+^<"^ CC1> (N)) + 0(0. 



This factorized form is already known from studies of the QCD hadro-production of heavy quarks 
(see also Ref. [38|]) and allows for a separate treatment of the resummation of threshold logarithms 
a" \n m jS (hard, m <2n) and the terms proportional to a"fi~ m In' jS (Coulomb, m <n). Note, that 
the matching constant g^\{pi s ) in the first case does not depend on the Mellin moment N, whereas 

in the second case gij{ C (& s ,N) does. In the following, we will focus on the computation of the 
one-loop hard matching coefficients which is the main new result of the present paper and which 



allows for the extraction of the expansion coefficients of gf?\ in Eq. (25 to NNLL accuracy. All 
explicit expressions are given in App. [B] 

Before doing so, we briefly like to comment on the resummation of Coulomb corrections, 
which accounts for the bound-state effects in the gluino pair and which exploits an effective 
description of QCD in the non-relativistic regime. To leading power, it is long known that the 
so-called Sommerfeld factor Ac sums the pure Coulomb corrections in momentum space (/3) 
corresponding to ladder diagrams |39) . One has jl4j , 
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71 OCc 



D, 



A c (x) 



(26) 



x 

exp(x) — 1 ' 

where we have introduced the quantity D\ = C\ /2 — Ca as a function of Ca = 3 and the quadratic 
Casimir operators Q of the final-state SU(3) co i or representation. For initial-state gluons these take 
the numerical values 



Q = {0,3,8} for I = {1,8,27}, 

Di = {-3,-3/2,1} for I = {1,8,27}, 



(27) 
(28) 



depending on the gluino pair being in the I = {1, 8, 27 {-representation. In the qq channel, only 
the anti-symmetric octet contributes at first order in the threshold expansion and we always set 
Ci = C A = 3, thus it follows Di = -3/2 from Eq. (28^ 



Ref. 31 

The notations for D\ vary in the literature: in Ref. 



(n)R« 



where n = I — 1. 



uses a different notation: The coefficients are denoted by D y HH , 

, d K l in Ref. |l6j and D Ra in Ref. |2lJ (where it is 



14 



explicitly given for top-quark production, which differs from the color configurations of gluino production) 
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A formal expansion of Eq. ( |26| in a s reproduces the NLO and NNLO pure Coulomb terms 
to leading power. However, the expansion does not converge due to the high inverse powers of 
j6 close to the threshold and those singular terms even cause a fixed-order expansion of Eq. ( [T9| 
beyond NNLO to be ill-defined. In the context of hadronic heavy-quark production this has mo- 



tivated detailed studies of the phenomenological effects of Coulomb resummation [40[. Methods 
and results for the combined resummation of threshold logarithms and the Coulomb corrections 
for heavy quarks have also been presented in |17 31 4lJ. The effect of Coulomb resummation for 
the total cross section is small, e.g. 0(1%) for the related case of heavy-quark hadro production. 

Let us now turn to the calculation of the necessary one-loop hard matching coefficients. To 
that end, recall that the NLO scaling functions fgg and /i» near threshold can be written in a 
factorized form with respect to the Born contributions as |l2| 



fi 
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where Cf s and define the hard one-loop constants to be determined from matching to a 
fixed order NLO calculation near the threshold. For the former quantity, we actually need the 
individual components Cf^, with respect to the final-state color configuration. Therefore, we 
decompose Eq. ( |29| as 



f (10) 



,(00) 
J gg,I 
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(31) 



where Di and C\ are defined in Eqs. (27 1 and ( [28] ). Note, that only the color-symmetric states 
1, 8, 27 contribute to the gluon-fusion channel, cf. Eq. whereas the anti-symmetric octet scal- 
ing function is suppressed in the threshold limit and, therefore, neglected. For quark-antiquark 
annihilation, only the anti-symmetric octet channel is considered, see Eq. d7j. Summation over all 
color configurations defines Cf in Eq. ( 29 as 



L i 



,'s'.'s' 



(32) 



with Ni given in Eq. ( |T3p . 

The analytic expressions for the one-loop matching constants Cf^ and = g can be 
extracted from Ref. 1 16 1, where the authors studied the QCD effects for a gluino bound state T 
with an invariant mass M. At the production threshold for gluino pairs, the differential cross 
section reads (cf. Eq. (5) in [16]) 

"''(SA'I 2 ,." 2 ,/';! = ^T / i(s,M 2 ,^,f/ 2 )^2lm{g I (0,M-2m ? + /r | }, (33) 



M 



dM 



where Tij^r.i denotes the hard scattering kernel, Q\ the Green's function of a non-relativistic 
Schrodinger equation, which accounts for the binding effects, and Y g the gluino decay width. 
Suppressing higher powers of f>, the hard function can be factorized as (cf. Eq. (19) in [16]) 



>T,I 



x-Born 
- r ij^T,l 
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where z = M 2 /s. Vy,i denotes the infrared-finite parts of the ultraviolet-regularized virtual 
corrections and TZijj the real corrections. In the threshold limit, these quantities are unaffected 
by the dynamics of the bound state formation. Thus, the explicit expressions can be taken over 
for the calculation of the gluino pair production cross section in the threshold region, where the 
(imaginary part of the) Green's function in Eq. ( [33] ) is set to one. The difference between the 
2 — > 1 and 2 — > 2 kinematics is encoded in the Born term J 7 ^ ^ t • In order to obtain the NLO 
hard kernels for gluino pair production, we simply have to replace the latter by our LO functions 



121 and ( 18 1. Setting the binding energy to zero, we further replace M by 2m<r, and thus z by p in 



Eq. (|34). 

In the related case of the QCD corrections to hadronic top-quark pair production, this proce- 
dure has been discussed in Ref. |42| , showing that the required NLO matching of the inclusive 
cross section to NNLL accuracy near threshold including the decomposition for color-singlet 
and color-octet states can be performed with the help of the NLO QCD corrections to hadro- 
production of quarkonium |43| (see also Refs. |44 45j). 

In the full MSSM, the number of both the quark and squark flavors, that enter the virtual NLO 
contributions, is given by tif = 6 and we set all squark masses equal. The one-loop matching co- 
efficients C! j depend on the chosen regularization scheme and we find in dimensional reduction 
DR adopted in Ref. & 
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(36) 



where we have defined L t ^ 



ln(m 2 /m~), as well as 

s 



'■) 



9(b 1 (r)-b i {r)+2b' 1 (r)) (1 - r) + 2 (9r - (Q + 1) ) c 5 (r 



(37) 



while Af 1 (r) is given in Eq. (31) of Ref. |16 . The results involve functions ai(r), bi(r), b\{r) and 
Cj(r) which have been defined as certain limits of scalar one-, two-, and three-point integrals in the 
Appendix of Ref. |16) . For the convenience of the reader, we give the explicit analytic expressions 
in App. [A] 

The quantity A^(r) in Eq. (36 i diverges in the limit r — > 1, but multiplication with the Born 



cross section gives a finite result for the NLO scaling function in the threshold approximation 
(which is actually zero). This is due to the factor (1 — r) 2 , which shows up in the expansion of 
Eq. (18 i. However, starting from 0(fi 5 ) (or 0(/3 3 ) in the color-summed result of Eq. |Toj)), there 
is no such factor and one would create an artificial divergence, if one inserts the full LO scaling 
function into Eq. ( |30) . On the other hand, it is clear that in the exact NLO result, the function 
Ag||(r) should also possess terms depending on /3. These cancel the problematic higher order 
terms when multiplied with the Born function. Moreover, they give additional contributions to 



the NLO cross section, which are not treated by the ansatz (34 1. Note, that in the gluon-fusion 
channel, Eq. (35' is free of artificial divergences and the function Af g (r) is well defined for all 
r > 0. 



Within the regularization scheme DR underlying Eqs. (35 i and (36 i, the strong coupling con- 
stant is understood to be evaluated at a hard scale where all squark flavors as well as the gluino 
contribute within the virtual corrections. Conventional QCD computations on the other hand 
employ the MS-scheme, and so does the program Prospino |l3|. In order to compare to the 
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numerical output of Prospino in the MS-scheme, we have to perform a scheme transformation 
and decouple the SUSY particles as well as the top-quark from the spectrum. The necessary 
change of the renormalization scheme for a s and the decoupling can easily be done with the help 
of formulae given in Ref. (46). Assuming that the top-quark is lighter than all sparticles, we have 
to add the following shift to the NLO scaling functions in Eq. (|5j 
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Jill 



- f {11) L 
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ln(r) + -L,, + 



n f 
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where n\ denotes the number of light (massless) quark flavors. This leads to the one-loop match 
ing constants in the MS-scheme with a total of nj quark flavors {nt = n\ + 1), 

r2\ H jr 
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C, (4 + ln(2) 
n f ( ln(2) - I 



71" 



+ 36-6 ln z (2) 



n 2 + %Af(r) r 



181 8 



ln 2 (2) 



43 
36 



TC 



3 L tg 



n 



f 
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(39) 
(40) 



We remark here, that the gluino bound state computation of Ref. [16] has been performed 
in the limit m t — > wherever possible. Thus, the dependence on m t in Eq. (35 1 for Cf 8 ^ R 



the DR-scheme and in Eq. (40 1 for C^ g MS in the MS-scheme is only logarithmic. The NLO QCD 
corrections to the inclusive cross section [12] coded in the program Prospino [13], on the other 
hand, account for the complete dependence on mt. 



m 



In order to cross check our analytic results in Eqs. ( p9[ ) and (40 1, we numerically extract the 
one-loop hard matching coefficients of the color-summed NLO scaling functions at their threshold 
from Prospino, cf. Eq. ( |32| . For our numerical analysis, we set the squark masses to 600 GeV 
and vary the gluino mass between 100 GeV and 2 TeV in steps of 100 GeV. The top-quark mass is 
set to 175 GeV in the on-shell scheme. According to Eq. ([9]), the LO scaling function of the gluon- 
fusion channel does not depend on any of the SUSY masses and the top-quark mass dependence 



of Cj^ MS has canceled in Eq. (39 1. Thus, for the case of gluon-fusion, we find agreement with our 



analytic result within a few per mill over the whole range of input values. For ^-annihilation on 
the other hand, we encounter a dependence on the mass ratio r in Eq. ( [TT] ) for equal squark masses 
and we expect deviations due to finite contributions proportional to the top-quark mass. We find 

differences between the expression for Cfi^ S based on Ref. Jl^jl in Eq. ( 40 1 and the result extracted 



from Prospino, which amount to the order of a few per cent especially for mass ratios r > 1. 
Altogether, this constitutes an important cross check, both of our derivation and of the original 
computation of the NLO corrections in Ref. |13"| . Moreover, as already noted, the gluon channel is 
dominant for collider physics predictions at the LHC. Therefore, we are able to provide extremely 
accurate predictions for the gluino pair production cross section in the threshold region. 

We are now in the position to present the NNLO cross section in the threshold limit exact to 
NNLL accuracy. All coefficients of the threshold logarithms In" at NNLO can be calculated 



from the resummation formula (20 1 with the exponent (21 1 after an inverse Mellin transformation 
and with the knowledge of the one-loop matching coefficients CJ l in Eqs. (39 1 and (40 1. Note, that 
at ) we only keep logarithmically enhanced terms proportional to powers of ln(p) as well as 
Coulomb corrections in the following. The two-loop matching coefficients j defined in analogy 
to Eqs. d29l and d30l are presently unknown and we set them to zero in the results for the NNLO 



cross section in the threshold limit below. The determination of the two-loop hard constants C| ] 
would require a complete NNLO calculation, which is beyond the scope of the present study. 
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For the gluon-fusion channel we obtain in this way the NNLO scaling functions in the thresh- 
old approximation as, 



r(20) 



,(00) 



(16/T 2 ) 2 

-8CS+ 3 



+ ^j-{~ 192ln2 ^ + ( 44 + 16Cl - f w ' - 192 M 2 )) MP) (4i) 

21n(2)^ +441n(2) +8Qln(2) -48 In 2 (2) - 16/r 2 j 
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+ 3841n(2) 
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262624 6584 



Ci + ( - 768 - 32Q + 1152 ln(2) )C\\ 



3 3 

ni (^p- + ^Ci - 1088 ln(2) - 32Q ln(2) + 576 In 2 (2) - 32n 2 ^ 

129696 ln(2) +5136Ciln(2) +96C 2 ln(2) - 133920 ln 2 (2) - 5184Qln 2 (2) 
62208 In 3 (2) + 5328 n 2 + 200Q n 2 - 72001n(2)7r 2 + 33264£ 3 - 48Q£ 3 



-16tz 2 D 1 (3 - 2Di(l + t; spin )) | rn(jS) + Cf, 



with Di and Q given in Eqs. (27 1, (28 1 and Cf^ T in Eq. J39} for MS(n, = 5). & denote the values of 
the Riemann zeta function. Likewise, for quark-antiquark annihilation, we find 



(20) 
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(167I 2 ) 2 

-297 + n, 
8192 



^ + ^ j 128 In 2 (/3) + (-138 + 4n, + 1281n(2)) ln(/3) + 12 
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+ 41n(2)) -1021n(2)+321n 2 (2) + 3 | r2 | 
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+256 In 2 (2) - -5-7T 2 ) + 38928 ln(2) - 35712 In 2 (2) + 12288 In 3 (2) 



13592 



9 



-7T 



9 

4480 



ln(2)7i 2 + + WdW 3 - 2D t (l + v spin )) \ ln(jS) + 
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-''11 



with the MS-scheme result for C^ 8 MS from Eq. (40 1. The results in Eqs. (j41| and (42 1 agree with 
Ref. 1 21 [ where the approximate NNLO cross section at threshold has been computed for massive 
colored particle production in an arbitrary SU(3) co i or representation of the final state. In partic- 
ular, they also contain subleading NNLO Coulomb terms and the non-relativistic kinetic-energy 



Eqs. |4l| 
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corrections, which do not follow directly from the resummed cross section | |20| , but have to be 
determined from matching to explicit NNLO computations p7|450| . The latter ones are given 
by terms proportional to D\(3 — 2D\(1 + P S pin)) ln(p') in Eqs. (41 1 and |42| and depend on spin 
configuration of the gg final state through the quantity v sp i n . For the gluino pair in a spin-singlet 
configuration as realized for the S-wave in the {1,8,27} symmetric color representations of the 
gluon-fusion channel it takes the value v sp i n = 0. For a spin-triplet as in the antisymmetric octet 
representation of the qq channel we have fspin — 

-2/3. 

For direct comparison, we also present here the one-loop matching coefficients Cfj and C^ 8 

the notation of Ref. (21) (cf. in Eq. (A.2) of that reference). Using Eqs. (29), Q, Q and 
SOI, we find for fif = ji r = ]i 



m 



Re[C 



W 1 

gg^gg i 
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2C™-6Ci-C A 32 - -7i 2 - 4 ln(2) L fl + L/ ( \-C l L ¥ , 



(43) 



Re[C 

L 1 



c ! c l^gg i 



2 C\\ a - 6 C A - C F I 32 
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71 



1 41n(2)L^ + Lj) + (y-|n / )L /( , (44) 



which displays an additional dependence on the renormalization scale due to the particular nor- 
malization of Ref. 1 21 [. Also note, that our choice n* = 6 corresponds to ft/ + 1 in the notation of 
Ref. (21). 



For completeness we briefly list all functions governing the scale dependence up to NNLO 
in the gluino pair production cross section. These can be computed by standard renormalization 
group methods (see e.g., |42|) in terms of coefficients f>\ of the QCD beta-function and the splitting 
functions P;y which govern the PDF evolution. For the hard functions f- 
have 



11) 



in the MS-scheme we 



11) 



Jij,\ 
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16 n 2 
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(00) 
kj,l 



,(00) 
Jik.l 



(45) 



where ® denotes the standard Mellin convolution, and repeated indices imply summation over 
admissible partons. The coefficients fi\ and P-p are taken in an expansion in powers of oc s / (4tt) 



as in the normalization (24 1 (see Refs. [32,331). The explicit expressions near threshold read 

^241n(8^ 2 ) -48-24 ln(2) 



r(U) 



r(00) 

hs.i 



16n 2 
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3 Hl 



(47) 



Recall, that ft/ denotes the number of light quark flavors. Likewise, at NNLO, the scale dependent 
part can be calculated by evaluating 



j n 
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where we have suppressed the color indices I. Inserting the threshold approximations of the 
splitting functions and the NLO scaling functions, we obtain the desired results for ij = gg,qq 
near threshold, which explicitly read 
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In Fig. |2j we plot the color-summed NLO and NNLO scaling functions. For comparison 
we also show the exact LO results given in Eqs. (|9jl and (10'. We use a gluino mass nig = 
750 GeV and squark masses = 600 GeV which correspond to r = 0.64. In the gluon-fusion 
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Figure 2: The scaling functions fgg and ffjq with i = 0, 1,2 and j < i in the MS-scheme. The masses are 
trig = 750 GeV and m~ q = 600 GeV.' 



channel, the dependence on r starts at NLO and is rather weak in the threshold region. For 
^-annihilation however, one has a stronger dependence already starting at LO. The NLO results 
and in Fig. 2 are exact. Similarly, as the results for NNLO scale dependent functions 

(21) (22) J 1 J 1 

fij and fjj are based on Eqs. (48' and (49 1, they are also exact at all energies even away 
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from threshold. For the genuine NNLO contributions j\- we plot our new results ( |4lJ | and 
( |42] |. The threshold approximation for the latter functions could, in principle, be improved by 
adding constraints imposed by the high-energy factorization, see (51) for related studies in top- 
quark hadro-production. However, given the large gluino masses currently considered, this is not 
immediately relevant for phenomenology at current and foreseeable LHC energies. 

It should be stressed further that the scaling functions in the gg channel exceed those of the 
qq channel by about one order of magnitude as shown in Fig. [2] Keeping in mind that also the 
parton luminosity at a proton-proton collider such as the LHC favors the channel with initial state 
gluons over the one with quarks in the TeV-regime, we conclude that gluon-fusion is by far the 
dominant source for ^-production at the LHC. Thus, the theory predictions of the inclusive gg 
hadro-production cross section are mainly governed by the gluino mass and are rather insensitive 
to the squark masses. 



4 Hadronic cross section 

Here we discuss the total hadronic cross section, which is obtained by convoluting the partonic 
scaling functions with the PDFs, see Eq. d3). For the numerical results we keep the threshold 
enhanced channels gg and qq at all orders up to NNLO, while we consider only the NLO contri- 
butions for gq, which are the leading contributions of this channel. As already discussed in Sec. |3j 
at the hadronic level, the gg-channel accounts for the largest part, whereas the contribution of the 
(^-channel is a few percent of the ^-channel, only. 

For reasons of convenience, the computation of the hadronic cross sections employs a grid 
in the nig - -plane for the scale independent scaling functions f^°\ ij = gg>(\^ig(\i which has 
been extracted from Prospino. This grid has already been applied in the numerical check of 
the one-loop matching constant c[ 8g ' and c[ q ^ (see the previous Sec. 3 1. For the hadronic cross 
section computation considered here, these scaling functions are used to calculate the exact scale 



dependent scaling functions f-^K ffj^ and from the renormalization group Eqs. ( jlsj , (48 

and (49b. Using these results together with the threshold approximation for from Eqs. (41 
and ( |42] defines the theory predictions at approximate NNLO (dubbed NNLO a pp r ox in the sequel). 

We work in the MS-scheme, which is implemented in Prospino with n; = 5 light quarks and 
an on-shell top quark with mass m t = 175 GeV fll3) . The masses of squarks and stops are set equal 
to the value = A/5nig so that the gluino is always the heavier particle. We use the PDF sets 
ABM11 NNLO (53) and MSTW2008 NNLO PDFs (52) irrespective of the order of perturbation 
theory. In Fig. |3} we present total hadronic cross sections for gluino pair production at the LHC 
for the cms energies 7TeV,8TeV and 14TeV at LO, NLO and NNLO. The width of the bands 
indicates the theoretical uncertainty due to a variation of the scale ]i in the range \nig <]i< 2m g. 
The increase in the predicted rates due to the approximate NNLO corrections of the order of 
(9(15 — 20)% at nominal scales is clearly visible and cross section numbers for selected gluino 
masses are given in Tabs. JTj — [3j Over the plotted range of nig, the cross sections in Fig. |3] are 
decreasing over more than four orders of magnitude. 

In Fig |4j we show as an example the scale dependence of the hadronic cross section for = 
750 GeV, niq = m- ti = m~ ti = 600 GeV for the LO, NLO, and approximated NNLO cross section. 
The cross section with its uncertainty at the LHC with 14 TeV cms is 1.43±°;§pb, 2.16+°;^ p b, and 
2.56lQ Q7pb at LO, NLO, and NNLO, respectively. One observes a strong decrease of the scale 
uncertainty. The K factors are K NLO = ^nlo/^lo = 1-46 and K NNLO = c^nnlo/^nlo = 1-13 and 
the point of minimal sensitivity, where the cross section adopts similar values for all orders is at 
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Figure 3: Total hadronic cross section and its theoretical uncertainty at the LHC for three different center 
of mass energies (7 TeV(upper figure), 8 TeV (central figure), and 14 TeV(lower figure)) at LO (blue bands), 
NLO (green bands), and NNLO a pp r0 x (purple lines) as a function of the gluino mass. The masses of the 
squarks and the stop are set to = 4/5m^. At NNLO a pp ro x/ the theoretical uncertainty has shrunk to a 
small band. 
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Figure 4: Scale dependence of the total hadronic cross section at the LHC with 14 TeV for trig = 750 GeV 
and m.q = 600 GeV. The vertical bars indicate the total scale variation in the range [wz^/2,2?7Zg], the vertical 
dashed gray line in the middle of the figure indicates the cross section at the nominal scale y. — wig. 
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m g 


a 


(IX)) [pb] 




cr(NLO)[pb] 




a (NNLO) [pb] 


[GeV] 


x - 1 

x 2 


x = 1 


x = 2 


x = \ x = 1 


x = 2 


x=j x = 1 x = 2 


MSTW 2008 NNLO 


300 


46.839 


33.335 


24.383 


57.703 


51.178 43.701 


62.941 


61.264 


58.908 


400 


7.230 


5.080 


3.673 


8.986 


7.919 


6.707 


9.814 


9.591 


9.207 


500 


1.475 


1.026 


0.735 


1.855 


1.624 


1.365 


2.029 


1.990 


1.906 


600 


0.359 


0.247 


0.176 


0.457 


0.398 


0.332 


0.501 


0.493 


0.471 


700 


0.098 


0.067 


0.047 


0.127 


0.110 


0.091 


0.140 


0.138 


0.131 


800 


0.029 


0.020 


0.014 


0.038 


0.033 


0.027 


0.042 


0.042 


0.040 


ABM NNLO 


300 


29.433 


21.365 


15.930 


36.863 


32.778 28.213 


40.176 


39.540 


38.548 


400 


4.012 


2.916 


2.176 


5.053 


4.493 


3.869 


5.507 


5.485 


5.388 


500 


0.739 


0.539 


0.403 


0.933 


0.831 


0.717 


1.017 


1.025 


1.015 


600 


0.165 


0.121 


0.090 


0.209 


0.187 


0.161 


0.228 


0.233 


0.232 


700 


0.042 


0.031 


0.023 


0.054 


0.048 


0.041 


0.058 


0.060 


0.060 


800 


0.012 


0.009 


0.007 


0.015 


0.014 


0.012 


0.017 


0.017 


0.017 



Table 1: Numerical values for the gluino pair-production cross section at LHC with yfs — 7TeV and the 
PDF sets MSTW 2008 NNLO (52), ABM11 NNLO (53). The QCD predictions are given at LO, NLO, and 
NNLO accuracy and for different gluino masses and scales x = jif nig. 

about ]i = 0.35ntg. 

Finally, we compare the total cross section for the two PDF sets MSTW2008 NNLO (52) and 
ABM11 NNLO (53) in Fig.[5j These PDF sets obtained in global fits differ significantly in the value 
of the strong coupling constant <x s and the shape of the gluon PDF at large parton momentum 
fraction x, e.g., oc s (M z ) = 0.1134 ± 0.0011 for ABM11 and a s (M z ) = 0.1171 ± 0.0014 for MSTW. 
These differences are marginally compatible and their origin has been discussed for instance 
in Ref. (53) . As a result, the cross sections calculated with the ABM11 set are of the order of 
O(30 — 50)% smaller over the whole range of gluino masses, see also Tabs. [2] and [3] As it stands, 
the differences in these non-perturbative parameter are the largest residual uncertainty in gg- 
cross section predictions with direct implications also for exclusion limits on nig and m-g reported 
by the LHC experiments. 

5 Conclusion and Summary 

We have studied the QCD corrections for gluino pair production at hadron colliders at NNLO in 
QCD. With the computation of the hard matching coefficients at NLO based on recent results for 
the production of gluino-bound states (16), we were able to derive all logarithmically enhanced 
terms near threshold at NNLO. Our results allow for the evaluation of the resummed gg cross 
section to NNLL accuracy or, alternatively, for predictions at approximate NNLO accuracy at 
fixed order of perturbation theory. We have chosen the latter approach to illustrate the impact of 
our new results on the apparent convergence and the scale stability of the hadronic cross sections 
at the LHC. In summary, we were able to promote the predictions for the gluino pair production 
cross section in the threshold region to the next level of accuracy, now putting it on par with 
squark-antisquark pair production. 

In advancing from NLO to approximate NNLO QCD predictions, we have found a significant 
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m g 


a 


(LO)[pb] 




a 


(NLO)[pb] 


(7(NNLO)[pb] 


[GeV] 


x 2 


x = 1 


x = 2 


x - i 

x 2 


x = 1 


x = 2 


x 2 


x = 1 


x = 2 




S4STW 2008 NNLO 


300 


75.431 


54.300 


40.118 


92.498 


82.484 


70.931 


100.577 


97.919 


94.424 


400 


12.576 


8.941 


6.532 


15.511 


13.760 


11.745 


16.887 


16.503 


15.894 


500 


2.768 


1.949 


1.411 


3.445 


3.039 


2.577 


3.755 


3.682 


3.541 


600 


0.727 


0.508 


0.365 


0.915 


0.802 


0.676 


0.999 


0.982 


0.943 


700 


0.216 


0.149 


0.107 


0.275 


0.239 


0.200 


0.300 


0.296 


0.283 


800 


0.070 


0.048 


0.034 


0.090 


0.078 


0.065 


0.099 


0.097 


0.093 


900 


0.024 


0.016 


0.012 


0.031 


0.027 


0.022 


0.035 


0.034 


0.033 


ABM11 NNLO 


300 


50.103 


36.595 


27.433 


62.375 


55.683 


48.137 


67.796 


66.555 


64.864 


400 


7.444 


5.435 


4.071 


9.320 


8.313 


7.182 


10.132 


10.047 


9.853 


500 


1.486 


1.087 


0.815 


1.868 


1.667 


1.441 


2.030 


2.034 


2.008 


600 


0.360 


0.264 


0.198 


0.453 


0.405 


0.351 


0.493 


0.499 


0.495 


700 


0.100 


0.073 


0.055 


0.126 


0.113 


0.098 


0.137 


0.140 


0.140 


800 


0.030 


0.022 


0.017 


0.038 


0.034 


0.030 


0.042 


0.043 


0.043 


900 


0.010 


0.007 


0.006 


0.013 


0.011 


0.010 


0.014 


0.014 


0.014 



Table 2: Same as Tab. [T]for the LHC with Vs = 8 TeV. 



increase in the rates, with K-factors of the order of 0(15 — 20)% depending, of course, on the 
chosen squark and gluino masses. The residual scale uncertainty on the other hand is generally 
small, of the order of a few percent only, showing good perturbative stability of the result. The 
largest uncertainty in the current predictions for gghadro-production is due to the necessary non- 
perturbative input, i.e., the value of a s (Mz) and the shape of the gluon PDF, where differences 
between the PDF sets ABM11 and MSTW amount to the order of 0(30 - 60)%. The impact of the 
latter differences on squark and gluino searches at the LHC is dramatic and the implications for 
any exclusion limits on squark and gluino masses has not been addressed so far in experimental 
analysis. 
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A Scalar n-point functions 

Here, we give explicit expressions for the one-, two-, and three-point integrals defined in Ref. [16| : 

a t (r) = r(l + 21n(2)-ln(r)), ©1) 

b x {r) = 2 + 21n(2) -rln(r) - (1 -r) In |1 -r\, ©2) 
b 2 {r) = 2 + 21n(2)-ln(r) + Re[ v / l^ln(-(l- v / l^)(l + v / T^)" 1 )], ©3) 
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a (NNLO) [pb] 
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x = 2 


X = 1 x = 1 


x = 2 
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MSTW 2008 NNLO 


500 


25.197 


18.566 


13.995 


30.549 


27.593 24.090 


32.894 


32.225 


31.340 


650 


5.027 


3.665 


2.736 


6.121 


5.506 


4.778 


6.599 


6.485 


6.303 


800 


1.270 


0.918 


0.680 


1.557 


1.395 


1.204 


1.681 


1.656 


1.609 


950 


0.375 


0.269 


0.198 
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0.414 


0.355 


0.501 


0.495 


0.480 


1100 


0.124 


0.088 


0.065 


0.154 


0.137 


0.117 


0.167 


0.165 


0.160 


1250 


0.044 


0.031 


0.023 


0.056 


0.049 


0.042 


0.060 


0.060 


0.058 


1400 


0.017 


0.012 


0.009 


0.021 


0.019 


0.016 


0.023 


0.023 


0.022 


ABM NNLO 


500 


17.216 


12.852 


9.815 


21.146 


19.120 16.772 


22.759 


22.453 


22.042 


650 


3.127 


2.332 


1.778 


3.860 


3.486 


3.055 


4.155 


4.131 


4.075 


800 


0.729 


0.544 


0.415 


0.903 


0.816 


0.715 


0.972 


0.974 


0.966 


950 


0.201 


0.150 


0.115 


0.250 


0.226 


0.198 


0.269 


0.272 


0.271 


1100 


0.063 


0.047 


0.036 


0.078 


0.070 


0.062 


0.084 


0.085 


0.085 


1250 


0.021 


0.016 


0.012 


0.027 


0.024 


0.021 


0.029 


0.029 


0.029 


1400 


0.008 


0.006 


0.004 


0.010 


0.009 


0.008 


0.010 


0.011 


0.011 



Table 3: Same as Tab. [T] for the LHC with y/s = 14TeV. 



h(r) 


= 1 +21n(2) +r(l -r) -1 ln(r), 


®4) 


h(r) 


= 2 + 2In(2)+rIn(r)-(l + r)In(l + r), 


©5) 


h(r) 


= 2 + 21n(2) -ln(r), 


©6) 


Mr) 


= 2 + 21n(2) + (r- 1 - 1) ln|l -r\, 


07) 


b[(r) 


= -1 -r In |(1 -r)/r\, 


®8) 


b' 2 (r) 


= \{l-r)- x +r\n(r)(l-r)- 2 {l + r)- 1 , 


© 9 ) 


cx(r) 


= i Re[ - Li 2 ((5 - r)(3 + 2^1 - r - r)- 1 ) + Li 2 ((l - r)(3 + 2y/l - 


-r-r)- 1 ) (|A]10) 




-Li 2 ((5 - r)(3 - 2Vl - r - r)" 1 ) + Li 2 ((l - r)(3 - 2y/l 


-r-r)- 1 ) 




+ Li 2 (4 (5 + (-2 + r)r)- 1 ) - Li 2 (4r (5 + (-2 + r)r)- 1 ) 






+ ln(l + (1 - r) 2 /4)(ln(-5 + r) - ln(-l + r)) 






- ln(l - 4r(5 + (-2 + r)r)- 1 ) In(r) ] , 




c 2 (r) 


= l Re [-Li 2 ((3-r)(l + r)- 1 )+Li 2 ((l-r)(l + r)- 1 ) 


©11) 




+ Ii2 ((_3 + r )(l +r )-i)_U 2 ((_l +r )(l + r )-l)] / 




cs(r) 


= -1 tt 2 + \ (Li 2 ((-1 + r)(l + r)- 1 ) - Li 2 ((l - r)(l + r)^)) , 


©12) 


c 4 (r) 


= -(l+r)- x (l -In(2)+ln(l + r) - (1 + r) _1 r ln(r)) , 


©13) 


cs{r) 


= ^Re[Li 2 (-l/r)-Li 2 (l/r)]. 


©14) 



Note that in the limit r —¥ 1, all c, (r) apart from C4(r) simplify to — n 2 /8. 



19 



B Explicit expressions for the resummation formula 

Here we give the process-dependent matching constants of the general resummation formula p0| 
in the MS(n; = 5)-scheme at NNLL accuracy. To that end, we find it convenient to introduce the 
parameter N = Nexp(7g) and rearrange the terms in Eq. (21 1 according to 



G tj/I (N) = ln(N) • g}A\) +gUA) + ««4i( A ) + • • • ' W 



U]VJ ' Ai],lVJ ' "6;;, 

with A = a s p>Q InN, so that there are no terms proportional to Euler's constant 7e contained in 

A 



the final results of the coefficients g| i(A). With these conventions, we find for the hard constant 



g° ggA = l + fl s |4Cf I -192 + Ci(-8-41n(2))+241n 2 (2) + 127r 2 | (§2) 
+fl 2 |l6Cf! (-48 + 2Q (-l + ln(2)) +721n(2) -481n 2 (2) +3/r 2 



-Qf - (-823 + 1465 ln(2)) -4000 ln z (2) +2496 ln J (2) + 126/r 2 

-3441n(2)7T 2 + (1296 + 481n(2))^ 3 ^ + C 2 (321n(2) - 641n 2 (2) + 4/r 

4 
~9 



ni(ci (92-1161n(2) +48 In 2 (2) -3tt 2 ) + 1744 - 2560 ln(2) 



-1776 In 2 (2) - 624 ln J (2) - 102tt 2 + 1081n(2)7r 2 - 336£ : 



32 

+ — ( -8207 + 12260 ln(2)) -101344 In 2 (2) +66784 In 3 (2) -23040 In 4 (2) 
+3292tt 2 - 79921n(2)7r 2 + 5664 In 2 (2) n 1 + 204tt 4 + (35728 - 49392 ln(2)) £ 3 
+L f( ^192Cfj (-1 + ln(2)) + Q (-472 + 6641n(2) - 576 In 2 (2) + 48 zr 2 ) 

--ni (-68 + 2Ci(-l + ln(2))+921n(2) -481n 2 (2) +3tt 2 ) -21616 
+31888 ln(2) - 26304 In 2 (2) + 10368 In 3 (2) + 1332tt 2 - 17761n(2)7r 2 + 8064£ 3 

+Lj ^-16 (-1 +ln(2))n/ -2568 + 2568 ln(2) - 11521n 2 (2) + 144/r 2 ^ | , 
in the case of gluon fusion. For the qq channel, we obtain 

4j, 8s = l + « s {4C 1 \ + ^(-82-91n(2)+81n 2 (2)+47T 2 )+L,(l4-^| @3) 
+fl 2 | - yC^ 8n (41 - 571n(2) + 321n 2 (2) - 2tt 2 ) + ^-n, ^9460 - 133721n(2) 

+8400 In 2 (2) - 2496 In 3 (2) - 489tt 2 + 432 ln(2) n 2 - 1344£^ - 29920 In 2 (2) 

8 2 / 

+ — (-83453 + 127247 ln(2)) + - ( 79424 In 3 (2) - 20480 In 4 (2) + 5167/r 2 

-10428 ln(2)7T 2 + 5248 In 2 (2) tt 2 + 168/r 4 ^ + ~ (1793 - 1849 ln(2)) £ 3 
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-Lp {^C\\ a (-1 + ln(2)) + ^ n, (409 - 5531n(2) + 2881n 2 (2) - 18/r 2 ) 



+ - (-8283 + 11819 ln(2) - 86401n 2 (2) + 23041n 3 (2) + 502/r 2 - 4081n(2)7i 2 
+1792^ M -L 2 ( y MOn/ (-1 + ln(2)) + 245 + 641n 2 (2) - 245 ln(2) - 8/T 



The one-loop matching constants C^ MS and C^ g MS are given in Eqs. (39 1 and (40 1 and the variables 
Q and D\ in Eqs. (j27jl and (28 1. Note that t he p resently unknown two-loop matching coefficients 
Cj i have been set to zero in Eqs. (B 2 1 and (B 3 1. 

The Coulomb corrections depend on the color configuration of the gluino pair, but are inde- 
pendent of the initial state. According to Eq. p5), they are collected into 



4f(N) = l-« s 4D I 7r 2 y|+fl 2 |-^D 2 7r 4 N + D I 7r 2 y|(ln(N)(-44+^) 04) 



124 



88m(2) + -n, 5-61n(2) - L„ 44 



For completeness, we also give the coefficients g| j n the convention of Eq. (Bit. Intro- 



ducing the abbreviations Lf r = \n(jAj / fij) and Lg r = ln(4m~/ we obtain 
g\ = A^Po 1 (l - 2 ln(l - 2A) + ln(l - 2A) A" 1 ) , 



= ^ Vfr ( 2A + " 2A ) + \ lnZ (! - 2A ) 

+ A^Pq 1 (2A L fr + ln(l - 2 A) L^j 



05) 
®6) 



+ A, (2) /3 2 ( -2A - ln(l - 2A) ) + /3 1 ± In(l - 2A) ( D, w + d£ , ) , 



gg. 



i 

2 



A 



1 1 -^V n ( 1 " 2A ) + 



1 -2A/ 



2(1 -2A) 



l+ln 2 (l-2A) 



\ +A + ln(l-2A) + 

2 2(1 -2A) 



1 



1 -2A 



1 + ln(l - 2A) 



A 



(i) 



1 -2A 



AL 



A 



1 - 2A V 2 



ln(l -2A) 



2(1 -2A) 



+ A' 2 ^o 1 (2A L /r + L fr (l 



1 -2A 



+ ^o 2 A (3) 



1 + ln(l -2A) 



A 



2(1 -2A) 



2 2(1 -2A) 



2 2(1 -2A) 



. 07) 
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